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SUi^IARY 

The nonlinear cva*vature ejspressions for a twisted rotor olade cr i lesz^ 
undergoing transverse bending in two planes, torsion, and extension are 
developed. The curvature expressions are obtained using simple geometric 
considerations, in contrast to other methods described in the literature. The 
expressions are first developed in a general manner using the geometrical 
nonlinear theory of elasticity. These general nonlinear expressions are then 
systematically reduced to four levels of approximation by imposing various 
simplifying assumptions , and in each of these levels the second-degree non- 
linear expressions are given. The assumptions are carefully stated and 
their implications with respect to the nonlinear theory of elasticity as 
applied to beams are pointed out. The transformation matrices between the 
deformed and undeformed blade-fixed coordinates, which are needed in the 
development of the curvature expressions, are also given for three of the 
levels of approximation. Tiie presexit curvatui’e expressions and transfcrrxaoi-x'. 
matrices are compared with corresponding expressions existing in the literature. 
These comparisons indicate some discrepancies with the present results in the 
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nonlinear terms* The reasons for these discrepancies are explained. Since 
both the nonlinear curvature expressions and the nonlinear transformation 
matrices are needed to develop the nonlinear aeroelastlc equations of motion 
of flexible rotor blades * the effect w stability of the discrepancies has 
yet to be assessed. As a by-product of this study the controversy regarding 
whether the uncoupled extenslonal frequency of a rotating beam Increues or 
decreases with Increasing rotational speed is resolA^d. 
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INTRODUCTIOH 


AeroeXaatio it ability aasociated with flap-lag-toz^im^axial coupling 
of flexible helicopter rotor blades involves both linear and nonlinear 
coupling anong the degrees of freedom. If the second-degree nonlinearities 
have to be considered, the governing eg,uations of motion must include all the 
second degree nonlinear terms. 7o develop these nonlinear ei^uatlons a need 
thus arises both for the second-degree nonlinear e^qtressions for the bending 
curvatures and torsional curvature (sometimes called ’’torsion” or ’’total 
twist”) and the second-degree nonlinear transformation matrix between the 
dafomed and undafornad blada-fixad coordlnatas. The curvaturas are also 
needed in nonlinear static and dynamic analyses of tubes transporting fluids , 
propellers, etc. Love (ref* 1), using Klrchhoff's kinetic analogy, in con- 
junction with direction cosines, developed the linear curvature expressions 
for an inextensible curved and twisted beam. Hunter (ref 2), using a vec- 
torial approach, extended Love’s work to the nonlinear case without invoking 
the inextensibility assumption. Another vector extension of Love’s work to 
the nonlinear case was given by Loll and Itote (ref. 3). References 1 to 3 
ail treat the problem which includes initial bending curvatures and initial 
twist, A special case of this problem which is receiving considerable atten- 
tion in the literature is that of a flexible hingeless rotor blade having only 
initial twist (sometimes called "pretwist” or "built-in twist”). Hodges 
(ref, 4 ) developed the nonlinear curvattire expressions for a blade with zero 
pretwist following the approach given by Kovozhilov (ref. 5), Nonlinear 
curvature expressions for .n elastic blade were also given in references 6 and 
7 where they were obtained by solving a differential equation for the trans- 
formation matrix relating the blade-fixed coordinates of the deformed and 
undeformed blade, T!ie nonlinear cxirvatures were also given in reference £ 
where they were obtained from simple geometric considerations in combination 
with Kirchhoff’s kinetic analogy. 


In t 


3 paper the .':.pproach empioyel in reference b will be extended to 


obtain tae 
of exastic 


renerai nonlinear curvature exp'ressions using the nonlinear theory' 
ity. In general, the nonlinearity of the equations of the theory c 
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e^.astlclty can haya both gaonatrical and physical origin. Qeomatrio non- 
linearity is associated vith the necessity to consider the defomed con- 
figuration to write the equilibrium equations and the need to include non- 
linear terms in the strain-displacement relations. Physical nonlinearity is 
associated with the necessity to consider the relations between the eoBq>onents 
of stress and strain as nonlinear. In the present development only geometri- 
cal nonlinearity is considered. The general large defomatlon ejgpressions 
for the curvatures developed herein are systematically reduced to four le^'fls 
of approximation by isqposlng various sisq>llfying assumptions (figure 1). The 
levels of approximation addzwssed are: (1) the particular case of large de- 

formations in which the elongations and shears are less than iinlty with no 
restrictions on the rotations i (2) the first case of small deformations in 
which the elongations and shears are negligible compared to unity with no 
restrictions on the rotations; (3) the second case of small deformations in 
which the elongations, shears, and rotations are negligible compared to unity; 
(4) the classical linear case of small deformations in which tne elongations, 
shears, and rotations ai’e negligible cot^iared to unity and the squares and 
products of the rotations are neglected compared to the strains. Here, and 
in the subsequent discussions, the terms “'elongations', "shears”, and rota- 
tions" have the same meaning as in reference The curvature expressions and 
the transformation matrices will be compared with corresponding expressions 
existing in the literature wherever possible. 

It should be remarked that , for convenience , the case of small deformations 
I is obtained as a special case of the general case of large deformations rather 
than from the partic\U.ar case of large deformations addressed herein. This con- 
sideration is also reflected in the form of the block diagram given in figure 1. 

In deriving curvature expressions for a deformed blade or beam the need 
arises to en^jloy Eulerian-type angles to effect a transformation between de- 
formed and undeformed blade- fixed coordinates. If nonlinear curvatures are 
being developed these angles must be treated as finite rotations. Since 
transformation matrices corresponding to finite angles of rotation are not 
commutative, the order in which the rotations are imposed is important, A 
preliminarj*- investigation of the nonlinear curvature esq^ressions as influenced 


k 


by tb« order in which the rotationeX transfoniatlona between the defomed end 
undeforaed blade- fixed coordinates are iapoeed was also iglL-tma in reference 8. 

Of the six rotational transforaation sequenMS possible reference 8 con- 
sidered twos fliqHlag-piteh and lae-flap-piteb. It was shown there that the 
torsional curvature eaqpression for an assusnd flap-lag-pitdi rotatimal sequence 
differed from the torsional curvature esqpression for a lag- flap-pitch rotm- 
tional sequence. Ihe present paper also examines more completely the effect 
of these two rotational trsnsfoznatlon sequences on the curvature esqpresslons. 
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S»<BOLg 


dinction coiinsi (1 ■ 1, 2, 3) 

coordinate along defomed elastic axis 

elastic deformations of arbitrary point- on elastic 
axis in radial, edgewise, and flapwise directions, 
respectively 

inertial axes in hul plane with origin at hub 
center-line 

blade-fixed axis system which translates with 

respect to x y 2 
000 

blade-fixed axis system, after deformation, which 

translates with respect to x y 2 

o' 0 0 

blade-fixed orthogonal axis system in defonned 
configuration obtained by rotating xyz; x^-eixis is 
tangent to the deformed elastic axis 

Eulerian-type rotational angles between xyz and 
* 3 ^ 3*3 

built-in twist angle (initial twist; pretwist) 

extensional component of Green’s strain tensor 

twist about deformea elastic axis 

torsional curvature (total rotation rate abcur x, 
axis) and bendin,! curvatures, respective ;; 


Spttcial notation: 

( f 

4arlTatii» with raaptct to 4( 

)/di 

( )• 

darivatlv* vlth raapact to x, d( 

)/d» 
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AltALmCAI. fiSVBLOPM£»t 


In thli ■teflon th« ftntrnl nonllnnnr eurmturt nxprntiiont and tha 
aasoeiitad tranafoimtlon aatrioat batwatn tlui dafonaad and undaforaad blada 
eoordinataa for i>. ballooptar rotor blada vi U firat ba darlvad for tba oaaa 
in vhlob tbara ara no raatrlotlona <» tha alongationa» ahaarst and rotations. 
Slagle gaoaatrlc oonaidaratl^a will ba uaad. Ihe genaral e3q>raaalons will 
than ba ayatamatioally raduoad to four lavala of approxiaaition. In aach 
of thesa lavala* tha sacond-dagraa curvature axpraaaions will be given. For 
tha three cases of small deformation oonsidarad* the transformation matrices 
will also be given. Finally* the geometric development will be supplemented 
by a ooo^lementary derivation based on a direction cosine approach. 


Geometric Approach 


General case of large deformations . - A schematic representation of the un- 
deformed and deformed blade geometries associated with both a flap-lag-pitch 
and a lag- flap-pi ten rotational transformation sequence ts shown in figures 
2 and 3, respectively. ’The translational elastic deformations experienced by 
an arbitrary point on the elastic axis of the blade are denoted by u, v, w. 

The coordinate axes XYZ are inertial axes situated at the root end of the 

(nonrotating) blade; x y z u.e axes fixed to the blade at an arbitrary point 

o o o 

on the elastic axis of the undeformed blade. Before deformation x y^z are 

0 0 0 

parallel to XYS. Deformations u, v, w, and displace the triad to xyz 

and rotate xyz to x^y^z^ where the axis x^ is tangent to the deformed elas- 
tic axis. The rotation of ti.e triad to its fiaal position x„y,z may be 

j j 

expressed in terms of the Sulerian-type angles^ 3» ajid 8 as shown in 
figure t for a flap-lag-pitch rotational transfonaation sequence and in 
figure 5 for a lag- flap-pitch rotational transformation sequence. two 


figures reflect the use of a 
,:r,et>mefci 7 of a l-lade ele/:cnt. 


on uescription tc describe the deforwed 


* *, ^ 'r ^ r* ^ * - a H t • i . <t v» 4 a , -i-i • , * r- 


•nan 


^Angies of this t^.-ie are 


further iiccuase 1 in references 8 ana j. 


H 


teufor wiXl ts^loyad to doflno tho •triina in t«jwHi of tho d*fonB*tloni» 
Bof«rono«s 2, 8, and 9 shovad that tha nonUnaar eurvatura axpraaaiona ara 
dapendant on tha ordar in vhieh tha rotational transfoxmationa hatvaan tha 
dafomad and imdafornad hlada-fixad eoordinataa avn tnpoaad« In tha pvaaant 
davalopnant tha handing and toraional eurraturaa eorraipondinf to tha two 
rotational tranafOraation aa^uanoaa, flap-lag»piteh and laft*f lap-pitch » ara 
darivad for tha ganaral eaaa of larga dafotaaelona uaing aiapla gaoaatrie 
conaidarationa and than apacialiaad to four Uvala of approxiaation, 

plap>lag-pitch transfomation aa<iuanca{ For thia rotational tranaforaa- 
tion aaquanca Cfigura h\ tha rotatlona ara ia^oaad aa follow: 

1« A poaitiva rotation 8 about tha nagativa y axia raaulting in 

2. A poaitiva rotation C about tha axia raaulting in 

3. A poaitiva rotation 0 about tha Xg axis resulting in 

Figure 1: reflects the use of Kirchhoff s xinatic analogy. This analogy (see, 
e.g, , refs. 1, 8, and 9) states that if the origin of the x,y_z_ coordinates 

Z ^ ^ 

moves along the deflected and twisted elastic axis with a unit linear velocity 
and the and axes occupy, at each instant, positions corresponding to 
the normal and binormal directions of the deformed elastic axis*, then the 
angular velocities of the system with respect to the xyz system are 

expressed in the same way as the corresponding bending and toi'sional curvatures. 
One then merely has to replace the time derivative in the expressions for the 
angular velocities by a space derivative with respect to the (curvilinear) 
coordinate s, along the elastic axis of the deformed blade. The space de- 
rivatives of the rotation angles with respect to s^, 8 , C , and 6 , as shown 
in figure U, assume the role of the angular velocities in application of the 
enel-ogy. The bending curvatures and , and the torsional curvature 

si’e obtained by projecting C"*, and O"*” along the axes and 


*?cr ol-e general case of large deformations, the principal axes of the deforiTied 
blade do not coincide- with the normal and binormal -U roctions of the ueformed 
elastic tCvis \,see refs 5 and 10). 


9 



by 




■ e'*' - ^ sin 5 

*3 
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»• rot.tlontl tr«.form.Uon Dctmen the i «,d *ye tried, (fisur. 2) eu 
be written in tenna of direction ooiinsa as"^ 
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For the flap-xag-pltch sequence being considered the explicit fora 


'2. is 


of equation 
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cot 6 cot C 


-tin c C 08 6 cot & 
-tin B tin $ 


tin i 

cot c cot 0 


cot i fin 0 


cot S tin d 
-tin t sin 0 cot 9 


■in r, cos S tin 9 
•tin 6 cot 8 


( 3 ) 


-cot 5 tin 8 


■in t sin 6 sin 8 
+cot $ cot 9 
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.4^ 
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the taek of expreeeing the rotation angles S, o and their 

*3 ♦’ “"1 ® Ji- 

of Xj vlth respect to xyz oar. ce vri-.ten 

: t’ei*ence as 


- 1/2 

' il+u') 
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REPRODUCIBILITY or Hit* 
I®IGINAL PAGH' IS PiKilt 



From tquatlons 


■in 3 ■ 


Ml 


V 


i ♦ 2c - V 


'.T 


•in C ■ 



cos 0 


♦ 2e - v'^ 



( 5 ) 


where c is the extensional component of Green's strain tensor and is given 

by 


e ■ u' + 1/2 (u’^ + v'^ + w'^) 


( 6 ) 


From th<^ defi j. ,'.on of Green's strain tensoi-, the derivatives with respect to 
are related to the derivatives with respect to x according to 

+ - 1/2 

( ) * (l+2e) ( ) ( 7 ) 

The third rotation angle 9 specifies the orientation of the y^^^ axes with 
respect to yg^o* ^is additional rotation is due to torsion of the blade, 

'i, about the x^ axis in the absence of pretwist. Common practice in the 
rotor blade j.lterature is to combine the pretwist with the elastic torsion. 

Using this expedient, the third rotation angle 9 in the present development is 
given by 

G = e ^ (£) 

p“ 
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As s r*«iat of ths dofonastions, th* lift* dx bscomsi an altnant of 

arc ds, of tba dafonnad alaatic axia* In viaw of thiif tha aaq?ressiona given 
in aiiuatlona 4 ara tha diraetion ooainaa of tba tangent (x^ nxia) to tha da- 
fomad alaatic axia. The coordina,ta axaa ara allgnad along tha tan- 

gent. normal, and binormal diw»ctiona of tha dafomad alaatic axia. Thaaa 
coordinate axaa and tha trigonomatric relatione given by aquationa 5 are ahoan 
in figure 4. 

Lag-flep-pitch tranafomation eequencat Por thia rotational tranaforma- 
tion aeauenca (figure 5> the rotatione are inq?oaad aa foUauet 


1. 

A positive rotation 

C 

about the 

% axis resulting in 

2. 

A positive rotation 

3 

about the negat.ve y^ axis resulting in 


*2^2*2 



Xg axis resulting in 

3. 

A positive rotation 

0 

about the 


Using figure 5 , the bending curvatures uiy^ and . and^the torsional 
cViTvature , are obtained by projecting 3**^, C » and 6 along the 

^3^3^ 3 ***** given by 


(ji * Q* + sin3 
*3 



cos$ sin6 - 3^ COS0 
Z* cos 3 COS0 + 3'*’ sinO 


(?) 


For the lag- flap-pitch sequence , the explicit form of eqiiation 2 is 



» 


cos B cos ^ 

-cos ? Bln S sin £ 
-sin Z cos 0 


-003 5 sin 6 cos 0 
+8 in C sin 6 


sin S cos 0 

cos ; cos 6 
-sir. 5 sin 6 sin 6 

-sin ; sin 0 cos 0 

-CCS t, sin 0 


sin 0 


COB 0 Bln 6 



COB 0 COB 0 


As before • the direction ooeinee of vith respect to xys can be 

written using ref^irence 5 as 



-1/2 


ij. 

(l+2e) 

(l-i-u 


-1/2 


“l“ 

(i+ae) 

V* 


-1/2 


“l * 

(l+2e) 

w* 


From eqxiatlons 11 


) ■ COS0 cosC 

■ sine cos3 (11) 

• sin$ 


sin 3 ■ 


^Lmm 

1 + 2e 


cos 3 


fl ♦ 2e 


1 


- w»^ 


1 + 2e 


sin C ■ 



V' 



♦ 2e - w''* 


cos C 




1 » u* 


1 + 2e 


— ;t 

- w’ 


( 12 ) 


Again, the third rotation angle is taken to be given by equation 8. Similar 
to the flap-lag-pitch case discussed above, equations 11 give the direction 
cosines of the tangent (x^ axis) to the deformed elastic axis. As before, 
the x^y^z^ coordinate axes are aligned along the tangent, normal, and bi- 
normal of the deformed elastic axis. Ihese coordinate axes and tr.e tri*o- 
nometrle relations given by equations 12 are showti in figure 5* 
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The exact Donllneu* curvature esqpreaelons and rotational tranafomation 
matrices given above » which are hi^ly nonlinear, nmy he approximated to any 
desired degree by applying the binomial theorem and substituting the trigo- 
nometric e:g;>anBion8 for 0, C» and 6* In the following sections, the general 
expressions will he reduced to four level's of approximation. In each of these 
levels <mly terms through second degvee in u» v, w, and ^ will he retained. 

The transformation matrices to second degree will he given for the three oases 
of small deformation considered. 

Large deformations . - For the particular case of large deformations con- 
sidered herein the second-degree nonlinear ejqjressions for the curvatures are 
obtained hy substituting equations 5i 6, Tt and S into equation* 1 for a flap- 
lag-pitch transformation sequence and equations 6, T» 8, and 12 into equations 
9 for a lag- flap-pitch transformation seqxience, expanding the resulting ex- 
pressions in a binomial series, and retaining terms up to second degree in u, 
V, w, and The resiilting curvature expressions assume the form 

» (1 - u* - |v»2 - e’ + 4.'{1-U»>- w” v' 

^*3 2 2 pt 

K, “ (“v" *♦• 2u'v” + u'V’ ) COS0 , + sin9 , 

/c pt pt 

+ (v" -2u’v" - u'V) sine^^ + cosO^^ (13) 

a {v'-2u‘ v'' - u"v’)cos6 . - ov'' sin9 . 

pt pt 

+ (;r '~2u'w'' - u'V) sir-t‘ ^ + tfrw’' cos9 , 


for 


*.ab-' 


ir’-pj-tCA o^3quence 



Wjj - (1-u* - I v*^ - I w»*> e* + ♦•<!-«•) + v'V* 
w ■ (-w*'+2u'w" + u'V ) co»e + v"sin6 . 

y 2 ^ 

+ (v*' - 2u*v’' - u*'w*) + 4> v" cos (lU) 

\ ■ (v"-2u*v" - u'V) cosSp^ - <J) v" sinQp^ 

+ (w"-2u*w" - u"w' ) sin6p^ + 4» w" co*®p^ 


for a lag-'flap>plteh saquanca* 

Small deformations I» - The majority of engineering materials are elastic 
only for small deformations characterized by elongations and shears which are 
negligible compared to unity. This implies that the strains are also negligible 
compared to unity. Since the rotations are arbitrary, the shears can be neg- 
lected in compeu'ison with the rotations in determining the direction cosines 
of the fibers of the beam in the deformed state. Tbis implies that plane 
sections remain plane (ref. 5, page kj). With the assumptions of this level of 
approximation the principal axes of the deformed blade differ little from the 
directions of the normal and blnormal to the deformed elastic axis. More de- 
tailed considerations on these aspects are given in references 5 and 10* in- 
voking the assumptions of this level the trigonometric relations of eqxxations 
5 and 12, after expanding in a binomial series and retaining terms up to 
second degree in u, v, and w, can be cast into the form 


sinB*w’ cos3 a 1 - “ (u'^ + w'^) 

1 2 

sine » v' cose s 1 - “ v' 

d 

and 

1 2 

sinpasw’ cos3 - 1 - ~ w' 

sinC'^.v' cose ^ ^ ^ 


(15) 


(16) 
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It should he observed that the second-degree ejq;>ression8 given in equation 1$ 
do not satisiy the trigonosietrie identity 

sin^0 ♦ cos^S *1 

and the expressions in equation l6 do not satisfy the identity 

sin^C + cos^C * 1 

2 2 
unless u' is negligible compared to xmity or, equivalently, u' is neg- 

2 2 

ligible compared to both v* and v* . This Implies that, under the assuap** 
tion that the elongations and shears are negligible compared to unity, the term 
u' in the extensional conponent of Green's strain tensor given by equation 6 
must be discarded in beam applications. This fact has not been pointed out in 
the literature as far as the authors know. It is interesting to note that this 
requirement is analogous to the inextensibility assumption invoked by Love 
(ref. 1) while developing the linear {first degree) curvature expressions for 
a bean. However, inextensibility has not been assumed in the present develop- 
ment as a consequence of this requirement as far as the second-degree nonlinear 
curvature expressions are concerned. Invoking the assumptions corresponding 

to the case of small deforiAatlons I and imposing the additional requirement that 

2 2 2 
u' is negligible conpared to both v’ and w’ , substituting equations 5* 

6, T» and 8 into equations 1 and equations 6, T, 8, and 12 into equations 9, 

expending the trigonometric expressions in a binomial series, and retaining 

terms up to second degree , the curvature expressions simplify to 


o; 


0 ’ -»• 3 ' - w'V 

nt 


CO 


“ -w" (cos0 ^ - -p sin- + v”(sin0 . + 1 cos6 , ) 

pt pt pt yz. 


( 17 ) 


= V i cos 


pt 


InO , ) + v' {3irJ- . + s cosv ^ j 

pu p u 
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for & flap-lag»piteh rotational sequance and 



0* + ♦ v'V* 

pt ^ 

-w"(cofep^- (p + v” (8ln0p^ + 4> co®0pt^ 

v"(co80p^ - 4> ain0p^) + w" («in0p^ ♦ 4> coaO^^) 


( 18 ) 


for a lag- flap-pitch rotational sequence. This is the level .if approximation 
usually employed for developing the nonlinear aeroelastic equations of motion 
of a flexible rotor blade. 

Small deformations II. - This case is obtained from the case of small 
deformations I by imposing the additional ass\jaq>tions that the rotations are 
negligible compared to unity and that the elongations and shears are much smaller 
than the rotations. If the first additional assumption is imposed the second- 
degree terms in equations 17 and 18 can be neglected compared to the firsu- 
degree terms and each of the equations reduces to the I'irst-degree (linear > 
expression 


0) 


X, 


9' + cjt 

pt 


to * -w''co30 , ' v"sin6 , 
7o P^ 


0 ) * v*’cos0 , + w'’sin8 . 

a, pt pt 


( 19 / 


It should be noted that the assumption that the elongations and shears are much 
smaller than the rotations has not been explicitly imposed. This assumption is 
required, however, in order to insure that the principal axis directions 
of a cross-section of the deformed blade still differ little from the 
directions of the normal and binormal of the dei’ormed elastic axis. It shell- 
also be pointed out that this is the level of approximation usually employe- ir. 
elastic stability (buckling) problems. 
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Classical linear theory. - To arrive at this case froa the case of small 
deformations II the additional assumption Is made that the terms of second 
degree In the angles of rotation are negligible conqpared to the corresponding 
strain terms. Herein t the terminology "classical linear theory" means that 
the geometric relations between strains and displacements are linear. Since 
there are no terms involving sq.uares of rotations In eqiuatlons 19 the curvatures 
for the classical linear theory are again given by equations 19. It should be 
emphasized, however, that although the c\urvature expressions corresponding to 
the last two levels of approximation are the same, the strain components are 
different (see ref. 5). 


Direction Cosine Approach 


For completeness it is appropriate to Indicate how one would extend the 
approach used by Love (ref. 1) for developing the linear curvature expressions 
to the general case of large deformations. The rotational transformation 
between the and xyz triads has already been given in equation 2. 

Following the procedure in Love , the curvatures written 

in terms of the direction cosines of equation 2 as 


-1/2 
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^ V 3 * “l®3 
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( 20 ) 


- 1/2 

'■he I'actor (l + 2s) in equations 20 accounts for the fact that the deriva- 

tives are taken with respect to the undeformed coordinate x. 

'ihe direction cosines in equations 20 depend on the order in whicii the 

rc.-!*tional transformation sequence between :c/z and is imposed. For the 

0 J o 

t f c r::uit i on sequences ocnsidered ^bove , the enn’icit foims ••qunhi c-ns 


nre riven by equation 
lut'- flap-pi ten sequence. 


for a flap-1 a j;-pi tea sequence .'js.’ equal I- r. 2 i for : 
Tae angles 0 , C, and 6 which appear in equatiui’S ?, 


le 



and 10 art dtfintd in t^uatlons 5* 6 » and 12 and thua tbt dlrtction eoainta 

art known in ttma of u» v» v, and 8^. Tht eurvaturt tkprtttlona for 
tht gtntral eaat of largt dtfomationa art obtaintd by o^atltuting tht di* 
rtetlon coainta igqprq;»rlatt to both tht flap-lag»pltch and lag-flap«>pitch 
tranaforoMtlon atqutncta Into tauation 20. Aa rtnarktd tarlitr thtat gtntral 
txprtaaiona art highly nonllntar and aty bt approxiaattd to any dtairtd dtgrtt 
by allying tht binomial thtortm. In thia atetion tht gtntral taprtaaiona viU 
again bt rtduotd to four Itvtla of approximation and only atcond-dtgrtt ttroa 
in u» V, V, and ^ will bt rttalntd. 

For tht particular caat of largt dtformationa considtrtd htrtin the 
atcond-dtgrtt nonllntar taq^rtaalona for tht ciurvaturta art obtaintd by uaing 
eauatlona 3» 5* 6 , and 8 in equation 20 for tht flap-lag-pitch atqutnct and uaing 
tquatlona €, 8 , 10 and 12 in equation 20 for the lag-flap-pitch atqiitnce , ex- 
panding the resulting expreaslona, and retaining terms up to second degree in 
u, V, V, and <p. The expressions thus obtained are identical with the second- 
degree c;xrvature expressions given previotisly in. equations 13 and lU. 

If the assimptlons corresponding to the case of small deformations I are 
imposed, the second-degree expressions for the direction cosines associated vi-h 
a flap-lag-pitch rotational transformation sequence , from equations 2 , 3 , 5 , c , 
and 8 , can be shown to be 
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and those associated with a lag-flap-pitch rotational transformation sequence, 
from equations 2, 6, 8, 10, and 12, to he 

* 1 - (u»^ + v'^ +w'^) 

i. - v'(sln6p^ + i oo^e^^) - «'(coee^ - « sine^^) 
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it sliould be 'ju^ervei: th.it the ‘lirection cosines tiveii in eauatic; 
do not nfitisT' the ortiiOt-on.-bl lt.y properr^y 
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2 

unless VI* ^ Is aasuned negligiblt coapsrsd to unity o* » eq^uivslontly t u 
is MsuDsd nsgllglble conqptred to both v** snd w*^. This it the stmt con- 
clusion rtacbtd terUtr on tht basis of trigononotric identities. Again, tbtst 
results reaffirm that if one assianes that the elongations and shears are 
negligible oonqpared to unity the additional requirement that u is negligible 
compared to both v'^ and must also be imposed. Substituting equations 

21 and 22 into equations 20 and invoking the above additional assunqition the 
second-degree nonlinear curvature expressions given in equations 17 and 18 are 
again obtained. Ihe corresponding rotational transformation matrices between 
the deformed and undeformed blade-fixed coordinates are given by 
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(23) 


for a flap-lag-pitch rotational transformation sequence and 
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for a lag- flap-pitch rotational transformation ae<iuence. The case of 
deformations II is obtained by imposing the additional restrictions that the 
rotations are negligible compared to unity and that the elongations and shears 
are much smaller than the rotations. In this case, the second-degree teres in 
equations 23 and 2U can be neglected compared to the first-degree terns and 
reduce to the single relation 
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on i rthogoncli ty relations tv :eCv.-;-i degree. 
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nonline&r t%xm it alto dietatod by tht atai»q»tlon that tht rotatioot art 
ntgliiihlt o<HB|^artd to unity. 

Again, tht curvaturta eorrtaponding to olattioal lintar thtory art ^ 
aasit at thoat eorrtaponding to tht caat of anall daforaationa XI. 

It vu indicattd tarlltr that tht ateond-dtgrtt eurvaturt tag;>rtaaiona in 
rtftrtnota 6 and 7 vtrt ohtaintd by aolving a difftrtntial aquation for tht 
tranaforoation natrix vhich rtlatta tht bladt-fixtd coordinatta of tht dtforatd 
and mdtfomtd bladt. If ont toqploya that proctdurt uaing tht t rant foaraat ion 
glv.n in equnticn. 23 2k .olv.. for th. nnknovn.. u*,. .^3. 
(O22 obtain tht aant curvatiirt txprtsaiona givtn in tquationa 17 and 

l8. 

The ateond-dtgrtt nonlinear eurvaturt exprtsaions for a twiattd rotor 
bladt or beam have been tstabliahed above in two vi^s for four Itvtla of approxi- 
nation. The sane expressions can also be developed if foreshortening due to 
bending (ref. 8) is explicitly included in the axial defonaation. In the next 
section the curvature expressions and transformation matrices developed herein 
vill be compared with corresponding results erlstlng in the literature. 

COMPARISONS AND DISCUSSION 

The only general development in the literatiire for the curvature express- 
ion. , as far ns the authors know, is that of Hunter (ref. 2). Since reference 
2 considers initial (built-in) curvatures before imposing the elastic defor- 
mations and the present development combines pretwist with elastic twist and 
considers them as the last rotation, a direct comparison of the present re- 
sults with those of reference 2 is not possible. A comparison can be made, 
however, if the initic.! curvatures in reference 2 and the pretwist in the 
present development are set to zero. Under these conditions the present 
£ 'cond-degree nonlinear expressions given by equations 13 for a flap-lag- pitch 
sequence and equations ll for a lag- flap-pitch sequence agree with the corres- 
ponding; expressions obtained from reference 2 after making the necesc-iry no- 
tatioual changes. 
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StvtrftX invcitii^tora h«vt aacond^dtcrM no&linaup eurmurc 

taq^ratalont Mitaaing tBall d«fonii«tioas. Bodftft (rtf* li) darivad tliatc ax- 

praaaiona for ztro pratwlat# ‘Ihotiilti not atatad axplioitly» tha laval of 

i|q;>roxiffi«tion «Bq;>lo3rad ia aqulvi^nt to ‘Uw oaaa of aiaall dafontjiticm I iuirain* 

2 

Howavar, rafaxanda 4 diaeaxdad u* tliroutft an a priori ordarlng aehaaa 
rathar thtti througli a rigozoua argiaaant. In addition to tliaaa aaata^ticma, 
rafaranca U naad a par^4ally nonlinaar diaplacanant fiald in u, v, v, and d* 

Tha uaa of thla diaplacdnant liald maana that tha rotatiixial tranafomation 
matrix hatvaan tha daformad and undeformad blada-fixed coordinataa vaa partially 
linaarized. If ona ia davaloping a linear aet of cvirvatura axpraaaions, a 
linear displacement field in u, v, v» and d &nd a linear resultant rotational 
transformation matrix are sufficient and the order in which the rotations are 
imposed is not in^ortant. However, if one is developing the second-degree non- 
linear curvature eaipressions care must he taken to insure retention of all 
second-degree terms in the displacement field and in the resultant transfor- 
mation matrix and the order in which the rotations are imposed is important. 
Since reference 4 did not address the rotational transformation seiiuence ex- 
plicitly and has employed an incomplete displacement field, a meaningful com- 
parison with the present results can not he made. However , it is interesting 
to note that the second-degree curvature expressions of reference 4 agree with 
those associated with the flap-lag-pitch sequence given herein after making 
the necessary notational changes. In view of the use of a partially nonlinear 
displacement field, this agreement must he regarded as fortuitous. 

?.efe->^ence3 6 and J have also aadressed the problem of deriving nonlinear 
curvature exj)ressions for a hingeless rotor blade. Ko indication is given in 
reference 6 as to the level of approximation' employed. The assumptions made 

in reference T are equivalent to those here .n corresponding to the aa^e f 

2 

small deformations I. Tae term u’ vas discarded in reference 7 cr. the haslr 


of an ordering scheme. Both these 
transforvition sequence uet^/eeu th'’^ 
coc relit': '-es . hoierenoe n a priori 
the pretwl-v being sere, deference 
.r.;rvt.' tire :,u> being equel to (u ^ + 

I 't 


references used a lag- flap-pitch rc 
deformed and the undelcmed llr.de- 
' -entl flea the titrsicinal curvatur*:' 

7 also a priori identified the tor 
■ ) ' . The definition of 7 herein 


tat ion 
f: xod 
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•«m M that ffivtn in rafamneet 6 and 7 « Tha raaulta of tha oltad rafarancaa 
diaagraa vith tha eorraapondlni ton ional ounratura for tha lat-flap-piteb 
■aquanoa darlrad baraln and glvan by ^ ^ ^ banding 

ourvaturaa of rafaranoa 1 ura» boiravar» In agraaMnt vith tha oorraaponding 
onaa haraln. Aa a eonaaquanoa of thaaa a priori idantifieatlona tha lag-fXap- 
pitoh rotational tranaforsiation natrix of nfaranoaa 6 and 7 do not agraa irlth 
tha eorraapondlng ona givan by aquation 2h harain. 

Stcond^dagraa nonlinaar eurvatura axpraaaiona have baan davalopad in 
rafarance 3 for application to dynanie analyaaa of tubaa conveying fluida* 

Tha lavtl of approxisuition anqloyad in rafarance 3 la equivalent to tha caaa of 
anukll deformationa I herein* Since the rotational tranafornation aequence vaa 
not apf^>'l;'led a cosqtariaon vaa made with both aeta of the preaent reaulta after 
Betting the initial curvaturaa and tviat of nference 3 and the pretviat in 
the preaent reaulta to zero. Thia conpariaon shoved that the toralonal curva- 
ture the bending curvature of reference 3 expressed in the 

present notation a^ree with the corresponding results herein. The bending 
curvature ia not in agreejnent with the present resxata. The detailed 

development leading to the curvature expresslona given in reference 3 were 
presented in reference 11. An examination reference 11 reveals that the 
disagreement with the present results is a consequence of the linearization of 
the resultant transformation matrix between the deformed and undeformed blade- 
fixed coordinates while developing the nonlinear curvature expressions. This 
linearization is not consistent vith the assumptions of small deformations I 
although It is justified under the assumptions of small deformations II. How- 
ever, as already stated the nonlinear terms in vue curvature expressions under 
the assumptions of small deformations II have no meaning and must be discarded. 

Under the assumptions of classical li^^;ar theory the curvatures are linear 
and are again given by equations 19. ^he linear curvatures can be obtained 
from the nonlinear curvatures «»r the case of small deformations I by discarding 
the squares and products of the efastic variables u, v, w, and and their 
derivatives. The linear cmu-v n** > orprescions '•/utained in this manner in 
references 2, 3, 4, T, and 3 ar<- hi agreement with the corresponding ones 
herein. 
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REPRODuciBiLrrY or lllE 

original page is P(K)h' 


During the present develo|«ient er. Itvterestlng observation has been made 
regarding the u’ term In the extensloaul coa;ponent of Green's strain tensor 
(eq. 6). These observations have a direct bearing on an existing controversial 
issue regarding tdiether the uncoupled extenslonal frequency of a rotating 
beam Increases or decreases with Increasing rotational speed. Beference 12* 
while studying only the extenslonal vibrations of a rotating beam, used the 
extenslonal component of Green’s strain tensor in the form 

e - o' + i u’* (26) 

2 

but was able to discard the terms associated with u' in the development by 
invoking the small strain assumption. This can be formally shown by rearrang- 
ing equation 26 into the form 

2 

1 + 2e « (1 u' ) ( 21 ) 


If the strain e in equation 27 is assumed negligible compared to unity then 
u' must also be assumed negligible compared to imity. Based on this assump- 
tion u' must be discarded in equation 26. With this assmption reference 12 
concluded that the extenslonal frequency of a rotating beam decreases with in- 
creasing rotational speed. It should be remarked that in the mere general case 

2 

including flapwise and edgewise bending the conclusion that u' must be dis- 
carded can not be established on the basis of the above reasoning. However, 
in the present development this result was established in the general case on 
the basis of considerations related to satisfying certain trigonometric and 
direction cosine identities. More recently, reference 13 contradicted the re- 
sult of I'eference 12 and concluded that the extenslonal frequency of a rotating 
beam increases with increasing rotational speed. In this i'eference the coupled 


flapwise and extensd 
using Jreen'r stnir: 


onal equations of motion of a rotating beam were derives 

2 

, tensor including the u' term under the small strain 


assumptic 


:he extensional equation of motion was obtained as a Epe>: 


case of the coupled equations , reference 13 was apparently unable to recognise 
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2 2 

that u* mat 1>a naglaettd coaq^ad to w* * at Idantlfiad hatain by tha 
trigonoMtrie tad diraetim eotlna Idantitiaa. Thata conaldaratioaa ravaal 

9 

that tha fraciuaney iBoraataa vith ineraaainf rotational tpaad If u* it ra> 

9 

tainad and daeraatat with Ineraating rotational tpaad if u* it diteardad* 

On tha haala of tha obaarvatlont made in tha praiant davalopoantt undar tha 

9 

attun^tion of tmall itraint, tha tarn u* mutt ha diteardad in tha axtaneional 
eonq;>onant of Graan*t ttrain tantor in beam applicationt. Thut, tha pratant 
authort u!‘a in agraamant vith tha ratult of rafaranea 12 that* nndar tha 
astumptlon of small strains , tha uneoiq^lad axtensional fraqxiency of a rotating 
beam decraases vith increasing rotational spaad* 

CONCLUDING RDMABKS 

A davelopmant was prasantad for tha nonlinear aaqtrasslons for tha banding 
and torsional curvatures of a twisted helicopter rotor blade or beam undergoing 
combined flapwlse bending, chordvise bending, torsion, and extension. The 
general nonlinear expressions developed are valid for large deformations. 

These general expressions were systematically reduced to four levels of approxi- 
mation and in each of these levels the second-degree nonlinear expressions were 
given. In specializing to the particular case of small deformations in which 
the elongations and shears (and hence strains) are assumed negligible conipared 
to unity with no restrictions on the rotations, trigonometric and direction 
cosine anomalies were identified. To remove these anomalies, it was shown that 
one must additionally impose the requirement that the square of the first 
derivative of the extensional deformation on the elastic axis be negligible 
compared to unity and the squares of the bending slopes on the elastic axis. 

This fact has not been pointed out in the literature as far as the authors 
know. Using this fact, a controversy existing in the literature regarding 
whether the uncoupled extensional frequency of a rotating beam inci'eases or 
decreases with increasing rotational speed has been resolved. Furthermore, 
the second-degree nonlinear curvature expressions ond the rotational trans- 
formation matrices between the naformed end •.nderoracd blade-fixed coordinates 
were compared with corresponding expressiorxs in the literature wherever 
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potsll}!** coBq?*riions isdioatad sevaral dlfcrapaneiei vith the present 
resiilts in the nooUneer tense. Bie reeeene for these discrepancies were ex- 
plalned« 
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CLASSICAL LINiAR THfORY 


Figure 1.- Levels of ipproximation within the georoetric nonlinear theory of 
elasticity which are addressed. 
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Figure 2.- Schematic representation of undelormed and cIefor«ed blade for flap^lag^ltdi totai:10Ral 

trans format i on sequence , 
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